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The extra energy, negative energy and annihilation of energy by the subwavelength conservative 
systems that have a wave nature of light or matter (quantum) objects are predicted. The creating 
and destroying of energy break the energy conservation law in any subwavelength physical system. 
The paradoxical phenomenon is demonstrated in the context of extraordinary transmission of light 
and matter through subwavelength apertures [T.W. Ebbesen et al., Nature (London) 391, 667 (1998) 
and E. Moreno et al., Phys. Rev. Lett. 95, 170406 (2005)]. 
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The energy conservation law is the most important 
of conservation laws in physics. Conservation of energy 
states that the total amount of energy in a isolated sys- 
tem remains constant. In other words, energy can be 
converted from one form to another, but it cannot be 
created or destroyed. The energy conservation law is a 
mathematical consequence of the shift symmetry of time; 
energy conservation is implied by the empirical fact that 
physical laws remain the same over time. The energy 
conservation affects all physical phenomena without ex- 
ceptions, for an example, the recently discovered extraor- 
dinary (enhanced) transmission of light through subwave- 
length apertures in a metal screen JJ, 0, H, US IS S HI- 
The transmission enhancement is a process that can in- 
clude the resonant excitation of surfaceplasmons 
Fabry-Perot-like intraslit modes [1, @, 0| , and evanescent 
electromagnetic waves at the metal surface Q. In the 
case of thin screens whose thickness are too small to sup- 
port the intraslit resonance, the extraordinary transmis- 
sion of light or matter (electrons) is caused by the res- 
onant excitation of surface waves @, d, 0> 0]- At the 
resonant conditions, the system redistributes the elec- 
tromagnetic energy around the screen, such that more 
energy is effectively transmitted compared to the energy 
impinging on the slit opening. The total energy of the 
system is conserved under the energy redistribution. In 
the present paper, we predict the extra energy, negative 
energy and annihilation of energy by an ensemble of light 
or matter beams produced by an array of subwavelength 
apertures. The creating and destroying of energy break 
the energy conservation law in any subwavelength physi- 
cal system. The phenomenon, in particular, is associated 
with the extraordinary transmission without assistance 
of the surface waves. 

Let us first investigate the transmission of light 
through a subwavelength structure, namely an array of 
parallel subwavelength slits. The array of M indepen- 
dent slits of width 2a and period A in a metallic screen 
of thickness b <C A is considered. The metal is assumed 
to be a perfect conductor. The screen placed in vacuum 
is illuminated by a normally incident TM-polarized wave 



with wavelength A = 2ttc/lu = 2tt /k. The magnetic field 
of the wave H(x,y,z,t) = U(x)exp(—i(kz + Lot))e y is 
assumed to be time harmonic and constant in the y di- 
rection. The energy balance, which determines the trans- 
mission coefficient of the slit array, is derived by calcu- 
lating the power of light beams in the far-field diffraction 
zone. The EM beams produced by each of the inde- 
pendent slits are computed by using the Neerhoff and 
Mur approach, which uses a Green's function formalism 
for rigorous numerical solution of Maxwell's equations 
for a single, isolated slit [TlT |. The transmission of 
the slit array is determined by calculating all the light 
power P(A) radiated into the far-field diffraction zone, 
x<e[— oo, oo] at the distance z>A from the screen. The 
total per-slit transmission coefficient, which represents 
the per-slit enhancement in transmission achieved by tak- 
ing a single, isolated slit and placing it in an M-slit array, 
is then found by using an equation Tm(A) = P{X)/MPi. 
Here, P is the total power of M beams produced by the 
array, and P\ is the power of a beam produced by the 
single slit. Figure 1 shows the transmission coefficient 
Tm(A), in the spectral region 500-2000 nm, calculated for 
the array parameters: a = 100 nm, A = 1800 nm, and 
b = 5 x 10~ 3 A max . The transmitted power was computed 
by integrating the total energy flux at the distance z — 1 
mm over the detector region of width Ax = 20 mm. The 
transmission spectra 1m(A) is shown for different values 
of M. We notice that the spectra Tm(A) is periodically 
modulated, as a function of wavelength, below and above 
a level defined by the transmission Ti(A) = 1 of one iso- 
lated slit. As M is increased from 2 to 10, the visibility of 
the modulation fringes increases approximately from 0.2 
to 0.7. The transmission Tu exhibits the Fabry-Perot 
like maxima around wavelengths A„ — A/n (n=l, 2, ...). 
The spectral peaks increase with increasing the number 
of slits and reach a saturation (TJ^ ax « 5) in amplitude 
by M = 300, at A « 1800 nm. The peak widths and the 
spectral shifts of the resonances from the Fabry-Perot 
wavelengths decrease with increasing the number M of 
slits. Figure 1 indicates that enhancement and suppres- 
sion in the transmission spectra are the natural proper- 
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FIG. 1: The per-slit transmission Tm(X) of an array of inde- 
pendent slits of the period A versus the wavelength for dif- 
ferent number M of slits. There are three Fabry-Perot like 
resonances at the wavelengths A n ~A/n, n=l, 2 and 3. 

ties of an ensemble of independent subwavelength slits in 
a thin (b <C A) screen. The spectral peaks are character- 
ized by asymmetric Fano-like profiles. Such modulations 
in the transmission spectra are known as Wood's anoma- 
lies. The minima and maxima correspond to Rayleigh 
anomalies and Fano resonances, respectively. The Wood 
anomalies in transmission spectra of optical gratings, a 
long standing problem in optics [l2j ]. follows naturally 
from interference properties of our model. The new point 
is a weak Wood's anomaly in a classical Young type two- 
slit system (M = 2). Figure 1 shows the extra energy 
(T > 1), negative energy (T < 1) and annihilation of 
energy (T < 1). The creating and destroying of energy 
break the energy conservation law in the system of M 
independent subwavelength beams (slits). 

To clarify the results of the computer code we have 
developed an analytical model, which yields simple for- 
mulas for the diffracted fields. For the fields diffracted 
by a narrow (2a <C A, b > 0) slit into the region \z\ > 2a, 
it can be shown that the Neerhoff and Mur model simpli- 
fies to an analytical one. For the magnetic H = (0, H y , 0) 
and electric E = (E x , 0,E Z ) fields we found: 

H y (x,z) = iaDF^{k[x 2 + z 2 ] 1/2 ), (1) 

E x (x,z) = -az[x 2 + z 2 ]-V 2 DFt{k[x 2 + z 2 } 1 / 2 ), (2) 
and 

E z (x,z) = ax[x 2 + z 2 \- x l 2 DF\(k[x 2 + z 2 } 1 / 2 ), (3) 
where 

D = 4fc- 1 [[exp(ifc6)(aA - fc)] 2 - (aA + k) 2 }' 1 (4) 
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FIG. 2: Electromagnetic fields in the far-field zone, (a) The 
fields Re(E x (x)) (A and D), Re{H y (x)) (B and E), and 
Re(WE z (x)) (C and F) calculated for M = 10 and A = 
1600 nm. The curves A, B, and C: rigorous model; curves 
D, E, and F: analytical model, (b) Re(E x (x)) for M=l: 
analytical model, (c) Re(E x (x)) for M—5: analytical model. 

and 

A = F^(ka) + ^[F {ka)F?(ka) + F 1 (ka)F^(ka)}. (5) 

Here, F^, Fq, Fq and F\ are the Hankel and Struve 
functions, respectively. The fields are spatially nonuni- 
form, in contrast to a common opinion that a subwave- 
length aperture diffracts light in all directions uniformly 
[13| . The fields produced by an array of M independent 
slits are given by E{x, z) — J2m=i E m {x + mA, z) and 
H (x, z) = X^m=i Hm{x + mA, z), where E m and H m are 
the fields of an m-th beam generated by the respective 
slit. As an example, Fig. 2(a) compares the far-field dis- 
tributions calculated by the analytical formulas (1-5) to 
that obtained in the rigorous model. We notice that the 
distributions are indistinguishable. The field power P, 
which determines the field energy, is found by integrating 
the energy flux | (S) t | = \E x H* + E* x H\(c/16tt). Thus, 
the analytical model describes accurately also the energy 
balance in the system of M independent subwavelength 
beams. The model does not only support results of our 
computer code, but presents an intuitively transparent 
explanation (physical mechanism) of the extra energy, 
negative energy, and annihilation of energy in terms of 
the constructive or destructive interference of the M in- 
dependent subwavelength beams produced by the multi- 
beam source. The creating and destroying of energy, 
which are associated with the extraordinary transmission 
without assistance of the surface waves, break the energy 
conservation law. Notice that the array-induced decrease 
of the central beam divergence (Figs. 2(b) and 2(c)) is 
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FIG. 3: The per-slit transmission coefficient T(\) versus 
wavelength for the Young type two-slit experiment [la ]. Solid 
curve: experiment; dashed curve: analytical model. Parame- 
ters: a = 100 nm, A = 4900 nm, and b = 210 nm. 



relevant to the beaming light 14[, and the diffraction- 
free light and matter beams [l5| . The amplitudes of the 



beams (evanescent fields) can rapidly decrease with in- 
creasing the distances x and z. However, due to the 
enhancement and beaming mechanisms (Figs. 1-4), the 
array produces a propagating wave with low divergence. 
Such a behavior is in agreement with the Huygens-Frcsnel 
principle, which considers a propagating wave as a super- 
position of secondary spherical waves. 

It is now important to understand the energy bal- 
ance in the two fundamental systems of wave optics, the 
single-slit and two-slit systems. The major features of 
the transmission through a single subwavelength slit are 
the intraslit resonances and the spectral shifts of the res- 
onances from the Fabry- Perot wavelengths Q . In agree- 
ment with the predictions [f| , the formula (4) shows that 
the transmission T = P/P = (a/fc)[Re(D)] 2 + [Im(D)] 2 
exhibits Fabry-Perot like maxima around wavelengths 
A n = 26/ 'n, where P is the power impinging on the 
slit opening. The enhancement and spectral shifts are 
explained by the wavelength dependent terms in the de- 
nominator of Eq. (4). The enhancement (T(Ai)«&/7ra 
[lit ) 1S m contrast to the attenuation predicted by the 
model 0. At the resonant conditions, the system re- 
distributes the electromagnetic energy in the intra-slit 
region and around the screen, such that more energy 
(T > 1) is effectively transmitted compared to the en- 
ergy impinging on the slit opening. The total energy of 
the system is conserved under the energy redistribution. 
This mechanism is different from those based on the cre- 
ating and destroying of energy by the multi-beam (multi- 
slit) system. The Young type two-slit configuration is 
characterized by a sinusoidal modulation of the trans- 
mission spectra T 2 (X) [3, [13] ■ The modulation period 
is inversely proportional to the slit separation A. The 
visibility V of the fringes is of order 0.2, independently 
on the slit separation. In our model, the transmission 
is given by T 2 ~ J[F^ (x^iF^ {x 2 )]* + F} (x 1 )*iF ( }(x 2 )]dx, 
where x\ = x and x 2 = x + A. The high-frequency 
modulations with the sideband- frequency / S (A) «/i(A) + 
/2(A, A)~l/A (Figs. 1 and 3) are produced like that in 
a classic heterodyne system by mixing two waves hav- 
ing different spatial frequencies, /1 and f 2 . Although 
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FIG. 4: The per-slit transmission Tm(A) versus wavelength 
for the different values of A and M: (A) A = 100 nm, M = 2; 
(B) A = 500 nm, M = 2; (C) A = 3000 nm, M = 2; (D) A = 
100 nm, M = 5; (E) A = 500 nm, M = 5; (F) A = 3000 nm, 
M = 5. Parameters: a — 100 nm and b = 10 nm. There are 
two enhancement regimes at A < A and A>A. 



our model ignores the plasmons, its prediction for the 
transmission (T^wl.l), the visibility (VwQ.l) of the 
fringes and the resonant wavelengths A„ ~ A/n compare 
well with the plasmon-assisted Young's type experiment 
[lit (Fig- 3). In the case of 6 > A/2, the resonances at 
A„ = A/n can be accompanied by the intraslit resonances 
at A„ = 2b/ n. One can easily demonstrate such behavior 
by using the analytical formulas (1-5). We considered the 
TM polarization because TE modes are cut off by a thick 
slit. In the case of a thin screen, TE modes propagate 
into slits so that magneto-polaritons develop. Because 
of the symmetry of Maxwell's equations the scattering 
intensity is formally identical with E and H swapping 
roles. The described mechanism is not the only con- 
tribution to enhanced transmission. There can be also 
enhancement due to the energy redistribution by surface 
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1 1 SI ] - The surface waves can 
couple the radiation phases of the slits, so that they get 
synchronized, and a collective emission can release the 
stored energy as an enhanced radiation. This kind of en- 
hancement is of different nature compared to our model, 
because the model does not require coupling between the 
beams. 

In order to gain physical insight into the energy bal- 
ance in the multi-wave (M>2) systems, we now consider 
the dependence of the transmission Tm(A) on the slit sep- 
aration A. We assume that the slits (beams) are inde- 
pendent also at A — > 0. According to the Van Citter- 
Zernike coherence theorem, a light source (even inco- 
herent) of radius r = M(a + A) produces a transver- 
sally coherent wave at the distance z>irRr/\ in the re- 
gion of radius R. Thus, in the case of a subwavelength 
light source (A « A), the collective emission of the en- 
semble of slits generates the coherent electric and mag- 
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netic fields, E = J2m=i E m exp(iip m )mM Eiexp(iip) and 
HKiMHiexp(itp). This means that the beams arrive at 
the detector with the same phases ip m (x) ~ <p(x) (see, 
also Ref. Consequently, the power (energy) of the 

emitted light scales with the square of the number of 
slits (beams), P « M 2 P±. Therefore, the transmission 
(Tm = P/MP±) grows linearly with the number of slits, 
Tm«M. For a given M, the function Tm(A) monoton- 
ically varies with A (Fig. 4). At the appropriate condi- 
tions, the transmission can reach the 1000-times enhance- 
ment (M = Az/TrR(a + A)). In the case of B > Xz/wr 
(A>A), the beams arrive at the detector with different 
phases ip m (x). Consequently, the power and transmission 
grow more slow with the number of beams (Figs. 1-4). 
Notice, that according to the energy conservation law one 
should expect that the energy (transmission) would re- 
main constant with changing the slit (beam) separation. 
The creating and destroying of energy in a wave field, 
which are associated with the extraordinary transmis- 
sion, break the conservation law. In the case of A >> A, 
our model is in agreement with the conservation law and 
theories of gratings [20]. Our consideration of the sub- 
wavelength gratings is similar in spirit to the dynamical 
diffraction models 
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22j . the Airy- like model _[23[ . and 
especially to a surface evanescent wave model [8] . 

One can easily demonstrate the breaking of energy con- 
servation in any subwavelength physical system by tak- 
ing into account the interference properties of Young's 
double-source system. At the risk of belaboring the 
obvious, we now describe the phenomenon. In the 
far-field diffraction zone, the radiation from two pin- 
holes of Young's setup is described by two spherical 
waves. The light intensity at the detector is I(r) = 
\(E/ri)exp(ikri+ipi) + (E/r 2 )exp(ikr 2 +tp 2 )\ 2 = I1+I2 + 
2{IiI 2 ) 1 / 2 cos{[kri + ipi] — [kr 2 + </? 2 ]). The corresponding 
energy is W = / j(h + h + 2{I 1 I 2 ) 1 / 2 cos([kr 1 + Vl ] - 
[kr 2 + (p 2 ])dxdy. Here, we use the units cAt/8ir = 1. 
In conventional Young's setup, which contains the pin- 
holes separated by the distance A >> A, the interference 
cross term (energy) vanishes. In accordance with the 
conservation law, the energy is W = J J(h + I 2 )dxdy = 
Wi + W 2 = 2W , where W\ = W 2 = W . In the case 
of Young's subwavelength system (A << A; r\ = r 2 for 
any coordinate x or y), the energy is W = W\ + W 2 + 
2 J r(/ 1 / 2 ) 1 / 2 cos(( / 9i — tp 2 )dxdy. The first-order correla- 
tion term is the positive or negative extra energy. At the 
phase condition ipi — ip 2 = it, the interference completely 
destroys (W = 0) the energy. The system creates energy 
(W = 4 Wo) in the case of ipi — ip 2 = 0. The same phase 
conditions provide the creating or destroying of energy by 
quantum two-source interference (for example, see formu- 
las 4. A. 1-4. A. 9 24]). The phenomenon depends neither 
on the nature (light or matter) of the waves (continuous 
waves or pulses) nor on material and shape of the sub- 
wavelength apertures (1-D and 2-D apertures or fibres). 
There is an evident resemblance between our model and 



a Dicke superradiance quantum model [25j of emission of 
an ensemble of atoms. A quantum reformulation of our 
model, which will be presented in our next paper, help 
us to understand better why a quantum entangled state 
is preserved on passage through a hole array [261 ]. 

It should be stressed that energy in conventional phys- 
ical systems can be converted from one form to another, 
but it cannot be created or destroyed. According to our 
model, energy may be created or destroyed by construc- 
tive or destructive interference of waves (beams) only 
at the extremely particular phase conditions. The in- 
terference completely destroys energy if waves interfere 
destructively in all points of a physical system. The in- 
terference creates energy if waves interfere only construc- 
tively. The experimental realization of such phase con- 
ditions is practically impossible in conventional physical 
systems. We showed that the waves generated by the 
point-like sources separated by the distance A < A (for 
example, subwavelength gratings) can satisfy the phase 
conditions in the far-field diffraction zone. 

This study was supported by the Hungarian Scientific 
Research Foundation (OTKA, Contract No T046811). 
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